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1. Introduction 

The radiative corrections for beta decay have been intensively investigated by a num- 
ber of authors, and the prime issue for such studies has been to deduce the value of the 
Cabbibo-Kobayashi-Maskawa (CKM) matrix element Vud from nuclear beta-decay data. 
An accurate value for Vud is important for testing the unitarity of the CKM matrix. The 
most precise values of Vud have been obtained from the accurate data of super-allowed 
0+ — > O"*" nuclear beta-decays [1]. Neutron beta-decay measurements provide an alter- 
native method of determining V^d, a method which does not depend on the accuracy of 
nuclear models. Neutron beta-decay experiments also provide the most precise determina- 
tion of the axial- vector coupling constant, g^, which plays an important role in hadronic 
weak-interaction reactions including many astrophysical processes. Theoretically, pion 
beta-decay can also be used for determining Vud- Unfortunately, however, the currently 
available experimental data on pion beta-decay are not accurate enough to allow us to 
take full advantage of this merit, see e.g. Cirigliano et al. [2]. 

To extract an accurate value of Vud from neutron decay data, the theoretical expression 
for the neutron decay rate including radiative corrections must be known with sufficient 
accuracy. The usual convention is to decompose radiative corrections of order a into two 
parts, the "outer" and the "inner" corrections [3, 4, 5] . The "outer" correction is a univer- 
sal function of the electron energy, independent of the details of the strong interactions. 
The "inner" correction stems from short-range terms and hadronic structure effects. This 
hadronic structure dependence (and additional nuclear structure dependence in the case 
of nuclear beta-decay) causes uncertainties in extracting fundamental quantities like Vud 
from experimental data ^. 

In this communication we present the first calculation of radiative corrections to neu- 
tron beta-decay based on a low-energy effective field theory (EFT). EFT provides symme- 
try constraints required by the underlying theory and a systematic expansion scheme for 
the evaluation of the hadron current. As suggested by Weinberg [7], low-energy hadronic 
physics can be described by an effective field theory of QCD known as "chiral perturba- 
tion theory" (xPT). The effective chiral Lagrangian, refiects the symmetries and the 
pattern of symmetry breaking of the underlying QCD. For massless quarks the QCD La- 
grangian is chirally symmetric, but chiral symmetry is spontaneously broken generating 
the pions as massless Goldstone bosons. Since the u and d quark masses are very small 
compared with the QCD scale Aqcd, and since the finite pion mass generated by the quark 
masses is small compared to a typical strong interaction scale, it is reasonable to treat 
the explicit chiral symmetry breaking terms as small perturbations. is expanded in 
powers of Q/A^ <^ 1 where Q denotes the typical four- momentum of the process in ques- 
tion or the pion mass, m^r, which represents the small explicit chiral symmetry breaking 
scale. The chiral scale, A^ ~ ^T^fir — 1 GeV (/^r = 92.4 MeV is the pion decay constant), 
is associated with the "high-energy" processes that have been integrated out in arriv- 
ing at Cy^ and with pion loops. The parameters appearing in >C^, called the low-energy 

new calculation of these radiative corrections, obtained with the standard model of electroweak 
interactions, has been reported in a recent preprint [6] . The results however seem to differ markedly from 
the classic calculations of Sirlin et al. [3, 4, 5]. 
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constants (LEC's), effectively subsume tlie liigli-energy pliysics that lias been integrated 
out. In principle, these LEC's could be determined from the underlying theory, but in 
practice the LEC's are determined phenomenologically from experimental data. Once the 
LEC's are determined from appropriate empirical data, then represents a complete La- 
grangian up to a specified chiral order. Furthermore, starting from C^, one can develop, 
for the amplitude of a given process, a well-defined perturbation scheme by organizing the 
relevant Feynman diagrams according to powers in Q/A^. If all the Feynman diagrams 
up to a given power, i/, in Q/A^ are taken into account, then the results depend only on 
the LEC's up to this order, with the contributions of higher order terms suppressed by 
an extra power of Q/A^. 

Over the past decade xPT has been successfully applied to many processes; for reviews, 
see, e.g., Refs. [8, 9]. Chiral Lagrangians including the photon field have been developed 
and applied to, e.g., pion-nucleon scattering, see Refs. [10, 11]. Our present calculation 
of the radiative corrections to neutron beta-decay is an EFT based on the spirit of the 
chiral Lagrangian approach. Thus we write down an effective Lagrangian, appropriate to 
neutron beta-decay, obeying chiral symmetry and involving a minimum set of LEC's and 
use the Lagrangian to estimate the relevant amplitudes to leading, next-to-leading, and 
next-to-next-to leading orders (LO, NLO, N^LO) in the Q/A^ expansion. In fact, since 
the typical energy transfer of the reaction is much smaller than the pion mass, the '^Q/A^ 
expansion" here has a special feature to be explained in the next section. 

The results of our EFT calculation confirm the expression for the model-independent 
universal function derived by Sirlin [3]. Furthermore, our calculation provides expressions 
for corrections of order a to the angular correlation coefficients in neutron beta-decay. 
We will show that the short-distance phenomena including the model-dependent hadronic 
radiative corrections can be condensed into two LEC's, one relevant to the Fermi constant 
Gp and the other to the axial coupling constant qa- The values of these LEC's need to 
be determined by experiments. In order to have crude order-of-magnitude estimates of 
our LEC's, we also compare our results with the "inner" radiative corrections obtained in 
the standard calculations. Furthermore, we shall argue that, provided the LEC's involved 
in our calculation are of a "natural" size, the neutron-decay rate and angular correlation 
coefficients calculated here are expected to have a precision better than 10^'^. 

2. Effective theory for neutron beta-decay 

Since neutron beta-decay is a low energy process, it is natural to use here heavy-baryon 
chiral perturbation theory (HBxPT), see, e.g., Refs. [8, 9]. In fact the appropriate ampli- 
tude, however without radiative corrections, can be obtained from HBxPT calculations 
of muon capture on a proton, fi + p —> n-\- u, which have been carried out including N^LO 
correction terms [12, 13, 14, 15]. Neutron beta-decay, however, has a feature not shared 
by muon capture, namely several different expansion scales. In particular, the maximum 
energy release, AM = — nip — rUe = 0.782 MeV, is very small compared to the pion 
mass and the nucleon mass mjv = {iTip + ninj/'i. Correspondingly, if we denote by 
Q the typical four-momentum transfer of the process, Q ~ AM is also very small. We 
therefore introduce here a particular "Q/A-^" expansion in which Q, unlike most HB^PT 
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calculations, only represents Q. The chiral symmetry breaking scale, itIt^/K-^ ~ 0.14, 
will be accounted for separately. The nucleon recoil terms are governed by the scale 
Q/iTLiq ~ 0.8 X 10~^, and they are NLO corrections to the LO expression. The scale 
Q/iTL]^ ^ Q/A-^ is numerically of the same magnitude as a/(27r) ~ 10~^, governing the 
radiative corrections, which are our primary interest {a is the fine structure constant). 
Therefore, for our present purposes, we consider the a/{27i) and Q/rriN corrections to be 
of the same order. 

The relevant effective Lagrangian, Cp, for the neutron decay process reads 



NN-y 



+ c 



evNN: 



where £21/7 is the lepton-photon Lagrangian, Cnn-j describes the heavy nucleon interacting 
with a photon, and CevNN gives the effective V — A interaction between the lepton and the 
heavy nucleon current. Since the pion mass is much heavier than the typical momentum 
scale of the reaction, Q <^ m^, we suppress the pion fields of the chiral Lagrangian, C^, 
and in jCp we have retained only the interactions between the heavy nucleon field, lepton 
current, and photons. Later in the text, we will discuss the role of the pions in the present 
calculation. Thus one has, through LO and NLO, 
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where F^^^ — d^A^ — dj^A^ and is the covariant derivative of QED. The ^a is the gauge 
parameter and we choose the Feynman gauge = 1. The is the velocity vector of the 
heavy-baryon formalism, which we take as v'^ = (1,0), and is the nucleon spin operator 

o 

25*^ = (0, a). The isovector magnetic moment in the NLO Lagrangian is l^y^ l^v = 4.706. 
The quantities ei, 62, ey and ca are defined as the LEC's of the theory. The LEC's ei and 
62 are the a-order corrections related to the wave-function normalization factors of the 
electron and proton, respectively. The LEC's ey and ca are the a-order corrections to the 
Fermi and Gamow- Teller amplitudes, where we have factored out the common coefficient 

{Gf Vud)/V2- Those LEC's arc used to absorb infinities coming from the virtual photon- 
loops and take into account short-range radiative effects. We remark that some of those 
LEC's contain contributions from, e.g., gis for the one nucleon sector without leptons[10] 
and XiS for the meson sector with leptons[16] in xPT ^. As is conventional, the parameters 

o o 

of the initial Lagrangian, e.g. the Fermi constant Gf and the axial coupling constant 9a-i 
^Unfortunately the connection between the LEC's ey and ca and the gi and Xi of Refs. [10, 16] is not 
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(a) ^ (b) (c) 




(g) (h) (i) 



Figure 1: Feynman diagrams for neutron beta-decay up to order a. In diagram (a), the 
four-fermion vertex can represent either the leading-order (LO) or next-to leading order 
(NLO) vertex, the latter being a l/ruN correction to the former. The crosses on the 
electron and nucleon lines in diagrams (c) and (e) are vertices involving the LEC's, ei 
and 62, respectively. The vertex of diagram (g) is given by the LEC's cy and ca- 



are taken as the coupling constants in the absence of radiative corrections and in the chiral 

o 

limit, = 0. Thus in particular, we assume that the Fermi constant, Gf^ Gp — 1.166 

xlO~^ GeV~^, as determined from muon-decay. As we discuss in the next paragraph, 
higher order hadronic corrections, i.e., pion- loops, renormalize these "bare" couplings to 

o 

their physical values in the absence of electromagnetic effects, e.g. 9a~^ 9a- Furthermore, 
radiative effects give rise to additional corrections to the coupling constants, Gp and 
qa which depend on the process being considered. These radiative corrections will be 
displayed explicitly in the present work. 

We calculate the Feynman diagrams shown in Fig. 1, where the vertices are determined 
by the Lagrangian, Cp^ given above. Several remarks are in order on the diagrams in Fig. 
1. Consider first diagram (a), which does not involve radiative corrections. Diagram (a) is 
a tree-diagram for the LO and NLO amplitudes. As regards the LO contribution, one may 
wonder why we do not consider here the pion-pole diagram (not shown). The pion-pole 
diagram, which is responsible for the induced pseudoscalar coupling, formally belongs to 
LO and hence would be included in normal circumstances. However, the extremely small 
momentum transfer involved in neutron beta-decay {Q -C m^) drastically suppresses 

straightforward. The Qi and Xi originate in Lagrangians involving only subsets of the degrees of freedom 
considered here and thus generate radiative corrections to only particular vertices in the diagrams for 
neutron beta decay. Their contribution can be absorbed in ey and ca, but ey.ci would also contain 
contributions from the LEC's of a yet-to-be-calculated Lagrangian involving the nucleon, lepton current, 
and photons simultaneously. 
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the pion-pole diagram contribution. Due to the presence of the pion propagator and a 
momentum of order Q at each vertex, the pion-pole diagram scales like {Q/rriT^Y ~ 3x 10~^ 
relative to the dominant LO terms. The accuracy of our present treatment does not 
warrant the inclusion of this tiny pion-pole contribution, and we will not consider it in 
the main body of our calculation. In the concluding section, however, we will briefly 
discuss the LO pion-pole term and its radiative corrections. Diagram (a) in Fig. 1 also 
includes the NLO vertex coming from the nucleon recoil terms oc Q/rriN featuring in 
Eq. (4). Since we are treating the Q/rrifq and a/{2-K) corrections as contributions of the 
same order, we will discuss these recoil terms later in the text; however, in evaluating 
radiative corrections, we need not consider the recoil terms since these corrections would 
be of higher order ~ a/(27r) x iiyQ / {'^^n) ~ 10~^. At order N^LO there occur two kinds 
of contributions. Higher order recoil corrections scale as {Q/mjqY ~ 10~^ and therefore 
can be neglected. The remaining N^LO terms (diagrams not shown) come from pion-loops 
and the corresponding hadronic LEC's which would appear in HB^PT Lagrangian at this 
order, see e.g. Refs. [8, 9]. The pion-loop diagrams which generate terms proportional 
to i.e., terms representing the hadronic vertex form factor effects, can be neglected, 
since their contributions are suppressed by a factor of (Q/A^)^ ~ 10~^ relative to the 
dominant LO terms. The remaining contributions of the pion-loops, which contain terms 
proportional to (m^/A^^)^, renormalize the bare quantities such as the "bare" axial vector 

coupling constant 9 a- These (m-^/A^)^ terms and the corresponding hadronic LEC's are 

absorbed into the renormalized gA so that to N^LO order, gA = 9a [1 + C^((^7r/Ax)^)], see 
e.g. Eq. (4.50) in Ref. [8] or Eq. (50) in Ref. [12]. Radiative corrections to the pion loop 
diagrams are suppressed by a scale {rriTr/A^)^ ~ 2 x 10^^ relative to the leading radiative 
corrections, and therefore their contributions can be ignored in the present calculation. 

The above discussion indicates that, to the accuracy in question, we need only consider 
radiative corrections of the following type. Of the contributions topologically represented 
by diagram (a), consider those involving the LO vertex and evaluate all possible radiative 
corrections applied to these LO diagrams. Diagrams (b), (d), (f) in Fig. 1 are one- 
photon loop corrections for the electron propagator, the nucleon propagator, and the 
four-point vertex function, respectively. Meanwhile, diagrams (c), (e) and (g) represent 
the contributions of the counter terms, the Ci, 62, cy and terms, in the Lagrangian. 
These LEC's remove the ultraviolet divergence arising from the loop diagrams (b), (d) 
and (f). As is well known, the infrared divergences contained in diagrams (b), (d), (f) 
should be canceled by the infrared divergences in the bremsstrahlung diagrams (h) and 
(i)^, and we have confirmed this cancellation explicitly. 

3. The correlation coefficients and the decay rate from EFT. 

A general expression for the differential neutron decay rate dP is well known [18] for a 

case AA liereiii only the neutron is polarized, and in which the nucleon recoil and radiative 

■^Recently these bremsstrahlung diagrams have been studied by Bernard et al. for radiative neutron 
beta decay, n^p + u + e + 'y, in EFT[17]. 
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corrections are ignored: 



dr 



(27r)^ 



+h-[Ap + Bp^ + D 



Pe X 



(5) 



Here and pe {Ey and p^j) are the electron (neutrino) energy and momentum, n is 
the neutron spin polarization vector, j3 = Pe/E^, and a,b, A, B, D are the correlation 
coefficients. If we calculate diagram (a) in Fig.l in the LO approximation, and if we neglect 
the nucleon recoil terms in the phase space factor, then our calculation reproduces Eq. 
(5), and furthermore we recover the standard lowest order expressions for the correlation 
coefficients as given in [18]: 



-2g\ + 2gA 



B 



29I + 2gA 



(6) 



where qa is the physical axial coupling constant. The coefficient b in Eq. (5), which reflects 
the presence of scalar and tensor weak couplings, vanishes in our LO calculation, since 
our Lagrangian only contains the standard vector and axial vector weak interaction. The 
parameter D in Eq. (5) is related to time-odd correlations and hence it also should vanish 
in the LO calculation since our Lagrangian is T invariant. However, "induced" D terms 
can appear at higher orders. For instance, interference between the weak magnetism and 
the radiative corrections would generate a D term of order 10~^[19]. 

As we proceed to include the higher order radiative diagrams generated by the La- 
grangian of Eq. (1), we encounter infinities coming from the photon-loop diagrams in 
Fig. 1. In order to eliminate these infinities, we need to introduce counter terms with the 
corresponding LEC's in our Lagrangian. We renormalize these LEC's in the usual effective 
field theoretical method based on the dimensional regularization of loop integrals [8, 9] . 
The finite LEC's renormalized at the scale are given by:^ 



1 , , 3 

ev,A - 2(^1 + 62) + 2 



- - 7e + ln(47r) + 1 



+ 3 In 



rriN 



(7) 



This renormalization is adequate to remove all the infinities associated with virtual pho- 
tons which we encounter in this calculation. The differential neutron decay rate including 
the radiative corrections and 1/mjv corrections is found to be: 



dV 



F{Z,Ee) \pe\ E, 



dEpdQp^dfl 



"Pi. 



(27r)5 mn[Ep + E, + E,{p.p,)] 



\M\ 



(8) 



where we have retained the relativistic expression for the phase factor, and 



|M|2 = mnrripEeE, (l + ^ e{^) + ^ ^i' ^ 



*The convention for the dimensional parameter e used here is: d = 4 — 2e. 
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xCo(i?e)(l + m {l + (l + 1^ Ci(Se)/3 -P. 

+ (l + ^ '^"0 ^^'^^'^ ^ C'3(£;e)^ • • /3 + [C^{E,) + C5(£;e)/3 • Mn ■ P.] ■ (9) 

The explanation of the quantities appearing in this expression will be given below. We 
remark that, in order to arrive at this factored form, we have freely exploited the fact 
that terms of order (a/27r)^, {a/2T^){Q /uin) and {Q/rriNY can be ignored to the order of 
accuracy of our concern. 

In Eq. (8) the Coulomb part of the radiative correction has been extracted as an over- 
all factor and incorporated into the usual Fermi function -F(Z, E^) ~ 1 + {a/2Ti)5^^"'^^^ = 
1 + an/fSj for Z = 1. In Eq. (9) the finite LEG, Cy, featuring in the factor (1 + ^ Cy) 
subsumes those short-range radiative corrections to the Fermi constant G^p which have 
been integrated out in arriving at our effective Lagrangian. This point will be further 
discussed in the final section. The axial coupling constant, qa, which has been renormal- 
izcd by pion loops, is multiplied by short-range radiative corrections involving the finite 
LEG as well as Cy. For convenience, and to simplify the results, we incorporate this 
radiative correction to qa into qa defined by 



9A 



9A 



l + ^(e2-e!;) 



(10) 



and this qa has been used in Eq. (9). Recall that qa corresponds to the physical value, 
with all short-range radiative corrections removed. 

In Eq. (9), d^^^ represents the model-independent radiative correction to Gp, which 
depends only on the kinematics of the electron, while S^^ gives the model-independent 

— * — * 

radiative corrections to the coefficients of the angular correlation terms, /3 • pi, and n • (3. 
The explicit expressions for 5^^) and 5^^^ are: 



8^^ 



6^^ 



3 In 

+4 



mj 2 (3 




(11) 



(12) 



Here 



Spence function defined by: 



(m^ — rri^ + mg)/2m„ is the maximum electron energy, and L{z) is the 



Jo t 



(13) 
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The factor Co(-E'e) contains the recoil corrections to the overall rate. It is given by 



Co(i?e) = 1 + 



mjv(i + m) 



max 
e 



l + ~gi) + 2^iv~9A{(^' + l)Ee) ,(14) 



where we have used E^ — E'^°'^ — + 0(l/mjv)- The other coefficients Ci{Ee) {i 
1, 2, . . . , 5) are given by 



Cl{Ee) 

C2{Ee) 



a|l + 



rriN 
1 

rriN 



{g\ + 2ixy~gA + 1) ^ (^^ + 1) \%iiy~gAE, - ^E^^-g^igA + /^y )] 



1 + 3^1 
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Ee) + 
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£;e(//y - 1) ^2 + '^9Al^V + l" 
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^{9A + fiv) 
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2mN9A 



(e: 



2^^(1 + 3^1) 



+ 



(^^+l)(l + 3^i) 



where a, A, B are given by Eq. (6) with the substitution gA — * gA- It is to be noted that 
Eq. (9) exhibits angular dependences that are missing in Eq. (5). These extra angular 
dependences arise from the NLO contributions that have been included in the l/mjv 
corrections (which leads to Eq. (9)) but ignored in the LO evaluation (which leads to 
Eq. (5)). It has been a common practice to approximate the overall kinematic factor in 
Eq. (8) by applying an expansion in 1/mjv- If convenient, one could use the following 
approximation: 



rUpEl 



{Ep + E, + E^P-p,) 
where we have used 



(E, 



max 



Eef 



1 + {3E, - - ■ p,) 



E^ ~ {E'^"'' - Ee) 



rriN 



E. 



1 + ^(1-/3 -p.) 



rriN 



(20) 



(21) 



The angular dependence appearing in Eq. (20) needs to be considered simultaneously with 
the angular dependences contained in Eq. (9). 

The model independent radiative correction 5^^^ in Eq. (11) agrees with that obtained 
by Sirlin [3], while 5^^^ in Eq. (12) also agrees with the result reported by Garcia and 
Maya [20]. Wc note that recoil corrections have also been calculated in the literature 
using the conventional methods. For instance, Wilkinson [21] evaluated corrections to the 
decay rate and the correlation coefficient A, and Bilen'kii et al. [22] computed corrections 
to the decay rate and the correlation coefficient a. Furthermore, Holstein [23] considered 
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recoil corrections to all the observables for general nuclear beta-decays. Our results for 
the recoil corrections agree with those found in these previous studies. 

4. Discussion and conclusions 

As mentioned in the introduction, a prime issue in the studies of neutron beta-decay 
is to deduce the precise value of Vud from the experimental data. Another issue is the 
extraction of the value of qa from the data. We shall discuss here the significance of our 
present calculation in connection with these two issues. 

To obtain the actual numerical values of Vud and qa we need to know the values of the 
LEC's, Cy and e^, pertaining to the lepton-current nucleon-current vertex. These LEC's 
parameterize short-distance physics not explicitly included in the effective Lagrangian, 
jCp, and they need to be determined empirically using appropriate observables. This 
is an important line of studies for the future. Here, instead, we discuss simple order- 
of-magnitude estimates of the LEG Cy, which is the most important LEG in neutron 
beta-decay. Based on the general estimation of a photon loop diagram, one may expect 
the natural scale for this parameter to be of the order of {a/2n)ey ~ 2 x 10~^, with 
Cy ~ ln(me/A^). To obtain another rough estimate of Cy we may compare our result 
for the neutron decay rate obtained from Eq. (8) with Eq. (6) of Marciano and Sirlin [5] . 
Thus we introduce the premise 

- 41n f^) + 31n (^) + In (^) +2C + A,, (22) 

where uiw^t^z are the masses of the W, Z bosons and itia is the axial mass scale. As is 
customary, we define the Fermi constant Gp oi muon decay by absorbing the factor 1 + 
(3Q;/47r)ln(m^/m2) into Gp [24]. The contribution \n{'m'\y /mz) in Eq. (22) is actually the 
difference between the contribution of the Z-box diagrams ^ in neutron beta-decay and the 
contribution of the Z-box diagrams in muon decay. In Eq. (22), the major contributions 
to the right-hand side originate from the short-range virtual photon corrections to the 
Fermi transition from the weak vector and axial-vector vertices. The former gives the 
contribution, 3 ln(mi4//mjv), and the latter \n{m\Y /itla)- The C in the expression is 
the long-range mo del- dependent correction coming from the axial-current and anomalous 
magnetic moments of the nucleon, and is proportional to {^s 9 a) where ixs is the isoscalar 
magnetic moment of the nucleon. A value of 2C = 1.77 was found in Ref. [5] . In an HBxPT 
calculation, however, we have verified that a correction estimated from the diagrams of 
C is of higher order oc a/(27r) {Q/itln)'^ and can be neglected (see section 2). Finally, 
the Ag term, which includes a short-range strong-interaction correction, is very small: 
Ag ~ —0.34 [5]. In this connection, it might be of interest to decompose, following 
Cirigliano et al[2], our ey into two parts: Cy = efP + ey. The efP term describes the 
universal short-distance physics of electroweak theory discussed by Sirlin [24], while the 
Cy term describes short-distance hadronic physics. It is possible that the Ag term is 
associated with the ey term. The above considerations lead to a rough estimate, ey ^ 20, 

^The Z-box diagrams here refer to diagrams like the one in Fig. 1 (f), with the photon replaced by 
the Z boson; see Fig. 3 in Ref. [24]. 
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i.e., [a/(27r)]ey ~ 4 x 10~^, which is of a natural size as discussed above. The above 
comparison also leads us to expect that the dominant contribution to Sy comes from the 
short-range electroweak corrections. 

The LEG enters only as a radiative correction to qa in Eq. (10), and therefore it may 
seem that there is no significant motivation to remove the radiative correction :^(e^ — e^) 
from Qa defined in Eq. (10) and deduce the values of qa- Indeed, if we limit ourselves to 
neutron beta-decay, all the observables can be expressed using qa without referring to qa- 
However, since radiative corrections are specific to individual processes, there should be 
cases wherein the removal of ^(e^ ~ ^v) from qa has physical consequences and hence 

does play a significant role. A possible example is the Goldberger-Treiman relation, 
Qa'^n = f^ng^nN-i where g^^j^i is the pion-nucleon coupling constant. To elaborate on this 
point, it is useful to illustrate processes which necessitate the introduction of the LEG, 
e^. To this end, we consider diagrams containing the exchange of a pion (pion-pole) plus 
a virtual photon. These diagrams involve three distinct one-particle-irreducible vertex 
functions. The first type is a nucleon-nucleon-lepton-lepton four-point vertex in which 
a virtual photon couples to both the nucleon and the leptonic currents. This class of 
diagrams requires a counter term involving ca associated with qa- The second type is a 
lepton-lepton-pion three-point vertex wherein a virtual photon only couples to the pion, 
the pion-lepton vertex or the lepton, and this vertex is related to the pion decay constant 
/tt. Some of the LEG's arising from this type of diagrams can be found in the chiral 
Lagrangian considered by Knecht et al. [16]. These LEG's are also related to the "inner" 
radiative corrections calculated for pion beta-decay, see e.g. [25]. The third type is a 
nucleon-nucleon-pion vertex in which a virtual photon only couples to the pion, the pion- 
nucleon vertex or the nucleon, and this vertex is related to q^^n. The corresponding LEG's 
are the QiS appearing in Miiller and Meifiner's work [10]. To our knowledge, however, 
no systematic HBxPT study of the Goldberger-Treiman relation including the radiative 
corrections associated with each of the vertices has been done so far. In fact the radiative 
correction really has not been fully studied yet in the standard approach. Instead it 
has usually been assumed that e;^ ~ Cy, which makes the radiative correction to qa small. 
Such radiative corrections could contribute to the evaluation of the Goldberger-Treiman 
discrepancy, but there is clearly not yet enough information to determine whether they 
turn out to be significant in comparison with the chiral symmetry breaking term. 

As discussed in section 2, we have not included in our work radiative corrections in- 
volving the NLO vertex or pion loop diagrams. The former should be suppressed at least 
by a factor of jiyQ / {'^^n) — 2 x 10^'^, and the latter by a factor of (m^/A^)^ ~ 2 x 10~^ 
relative to the leading radiative corrections. Also omitted from our work are the isospin 
breaking effects, which are naturally incorporated in the N^LO heavy-baryon chiral La- 
grangian [8] not explicitly written in this paper. Recently, Kaiser [26] studied isospin 
violation corrections to CpVud using HBxPT and found that the isospin breaking cor- 
rections are of the order of 10^^. To the accuracy of our present concern, we can safely 
neglect the isospin violation corrections. 

We now summarize. Using the effective field theory for neutron beta-decay, we have 
calculated the decay rate of the neutron and the angular correlation coefficients including 
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recoil corrections and radiative corrections of order a. We have included all non-radiative 
terms through N^LO except those which are neghgible because of the extremely small 
value of Q for neutron beta-decay. Our results reproduce the model-independent radiative 
corrections and recoil corrections in the literature. The short-range radiative corrections 
of the earlier calculations are replaced in our theory by the two finite radiative LEC's, 
Cy and e% where ey affects Gp and the difference, ey — e% affects qa- Via comparison 
with the results of the existing model calculations, we have argued that the value of Sy 
is of a natural scale. An advantage of our EFT approach is the possibility of evaluating 
higher order corrections in a systematic way, and the possibility to parameterize the 
strong interaction dependent contributions in terms of well-defined LEC's, which can in 
principle be obtained from independent experiments. The next order corrections in the 
EFT for neutron beta-decay are estimated to be of the order 10~^ or smaller. Therefore, 
to the extent that the LEC's involved in the present calculation are of a "natural" size 
(as discussed above), we expect our expressions for the rate and the angular correlation 
coefficients to be accurate to better than 10~^. 
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